Dense granular flows: 
interpolating between grain inertia and fluid viscosity based constitutive laws 
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A scalar constitutive law was recently obtained for dense granular flows from a two-grain argument, 
both in the inertial regime (grain inertia) and in the viscous regime. As the resulting law is not 
exactly the same in both regimes, we here provide an expression for the crossover between both 
regimes. 
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I. INTRODUCTION 

The deformations of granular materials are usually cat- 
egorized into quasistatic deformations, dense flows and 
collisional flows p], 0] . 

In the quasistatic and dense regimes, the shear-rate 
dependence can be expressed 0] in terms of a parameter 
/ = T7, where T is the typical time for a single grain, ac- 
celerated by the pressure P (or force Pd 2 ), to move over 
a distance comparable to its own size d. For dry grains of 
mass m, the surrouding fluid can often be neglected and 
the pressure is resisted only by grain inertia ( "gi" ) @ : 
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i.e., 




(1) 



For granular pastes, if the grain inertia is negligible as 
compared to the viscosity of the surrounding fluid ("fv"), 
the force Pd 2 exerced by the pressure is resisted by a 
typical Stokes viscous drag force felt by a single grain 
moving at the velocity d/Tf v 0]: 
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i.e. 



jT fv = 7 



(2) 



T fv P 

Empirically, the dimensionless parameter / has been used 
in the form of a frictional constitutive law, in terms of 
the ratio between the shear stress r and the pressure P: 
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P 



(3) 



Empirically again, a universal function fill) seems to ac- 
count for existing data in both regimes 
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Recently, a two-grain argument (see Fig. [TJ, based on 
the same physics, led to a clear distinction between inter- 
grain approach and separation times [f| : 
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As a result, two different functions arose for the regime 
where the grain inertia dominates and for the regime gov- 
erned by the fluid viscosity: 



T + P 



T P 
P~~ 



(G) 
(7) 



They do not include the saturation (/i(/) — > [ii at large 
/) that is present in Eq. ^ and which reflects the onset 
of the collisional regime. 

In the present note, we derive an interpolation between 
both regimes represented by Eqs. © and ([7]). 



II. CROSSOVER FROM GRAIN INERTIA TO 
FLUID VISCOSITY 

The regimes discussed here are dominated cither by 
the grain inertia or by the fluid viscosity. In other or in 
future experiments, it may happen that the system lie 
in the crossover between regime "gi" and regime "fv". 
In the present section, wc discuss how it is possible to 
interpolate between both behaviours. 

Cassar et al. Q proposed a universal behaviour, given 
by Eq. (|4]), which they showed to be compatible both 
with the data in the "gi" regime and with that in the 
"fv" regime. Nevertheless, the dimensionless parameter 
/ in Eq. (fj| does not have the same meaning, as it is given 
either by Eq. ([1]) or by Eq. ([2]). Hence, in the crossover 
region, there is no obvious interpolation between both 
definitions of /. 

In the present approach, we not only have this diffi- 
culty with /, but we additionally have two different con- 
stitutive relations, namely Eqs. © and ([7]). We therefore 
need to go back to the equation introduced by Courrech 
du Pont for freely falling grains in avalanches, and 
used by Cassar et al. [H with a pressure P. Omitting all 
numerical coefficients: 
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FIG. 1: Schematic evolution within the granular material dur- 
ing shear, (a) One grain is transported quasistatically from 
position 1 to position 2, then falls into position 3 due to the 
applied pressure P. In a recent work 0|, we rather con- 
sider a pair of grains during the period of time when they 
are close neighbours (bl-3). First, the deviatoric (typically, 
shear) component of the stress, r, helps the pressure P es- 
tablish the contact (bl). Then, the contact rotates due to 
the overall material deformation (b2). Finally, the deviatoric 
stress overcomes the pressure to break the contact (b3). Be- 
cause the pressure is compressive in a non-cohesive granular 
material, the typical magnitude of the force transmitted be- 
tween both grains is stronger when the contact forms than 
when it breaks. 




FIG. 2: Raw predictions of the two-particle argument [f| in 
the regimes dominated by grain inertia and by fluid viscosity, 
see Eqs. flSE}- 



where v is the grain velocity, P d 2 is the typical force 
resulting from the pressure, and drjv is the Stokes drag 
force of the grain in the fluid. With vanishing initial 
velocity v(0) and position x(0), we can derive the grain 
position: 
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(9) 
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In the present two-grain approach, we will use this re- 



FIG. 3: Predicted response for a granular material in the 
crossover between the grain inertia regime and the fluid vis- 
cosity regime. r/P is the effective friction. 7 g i and 7f v are 
given by Eq. (|16[) . The "gi" and "fv" limits (r/P as a func- 
tion of J g i or as a function of Jf v ), represented on Fig. [21 are 
(approximately) given by the edges of the meshed surface. 



suit with a stress a = r ± P instead of P. Moreover, let 
us define: 



P* = 
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(10) 



(11) 



which are the values taken respectively by T and P at 
the crossover between both regimes, see Eqs. ([T]) and (0). 

With these notations, Eq. ([9]) provides the time at 
which a grain submitted to a stress a has traveled a dis- 
tance d to meet another grain: 
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(12) 
(13) 



The lifetime of a contact, given by Eq. can be 
rewritten using Eq. (d^J): 
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where n = rj P and where g is an approximation (precise 
up to within two percent) for the inverse of function / 
defined by Eq. ([13]): 
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(15) 



With the notations of Eqs. (|10Hlip . both limits of the 
dimcnsionlcss parameter / can be expressed as: 



I Ki = 



P/P* 



(16) 
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Hence, in order to show how it is possible to interpo- 
late between Eqs. ^ and (O, let us use Eq. (fT3| to 
plot /i = t/P as a function of both 7 T* / J P '/ P* and 
yT*/(P/P*). This is shown on Fig. [3 



III. CONCLUSION 

A recent two-grain argument provided two distinct 
constitutive laws for the rheology of granular media in 
dense regimes: one for the regime (gi) where the grain 
inertia is dominant, and one for the regime (fv) governed 
by the fluid viscosity. In the present note, we followed 



the same arguments to derive an interpolation between 
both laws. 

Further studies are now needed in order to (i) test 
whether the distinct predictions are compatible with the 
experiments, and (ii) develop theoretical arguments that 
could rationalize the saturation of the apparent frictional 
response as the collisional regime is approached. 
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